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Abstract 

When heavy nuclei collide, a quark-gluon plasma is formed. The plasma is subject 
to strong electric field due to the charge of the colliding nuclei. The electric field can 
influence the behavior of the quark-gluon plasma. In particular, we might observe an 
increased number of quarks moving in the direction of that field, as we do in the standard 
electron-ion plasma. In this paper we show that this phenomenon, called the runaway 
quarks, does not exist. 
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The phenomenon of the runaway electrons in the electron-ion plasma has been dis- 
covered by Dreicer QTJ] in 1959 and was later described by Gurevich [0 in 1960. They 
were interested in what happens to the distribution function of the electrons in plasma 
if there is a weak external homogeneous electric field applied. The somewhat unexpected 
discovery was that regardless of how weak this field was, the distribution function of the 
fast electrons was found to be strongly distorted. 

The essence of the phenomenon lies in the fact that the fast electrons travel large 
distances on the average between collisions. While the slow electrons will give away all 
the momentum they acquired due to the motion in the direction of the electric field at the 
instant of collision, the fast electrons acquire a lot of energy which can not be given away 
by collisions. The faster they travel, the less important the collisions become for them. As 
a result, we observe the flux of fast electrons flowing in the direction of the electric field. 
The bulk of the electrons which are still distributed according to Maxwell distribution 
play the role of the "heat bath" with the spring of the electrons in the direction of bigger 
energies in the momentum space. 

Recently it has been remarked by I.Kogan [Q that somewhat similar physical setting 
is achieved in a completely different situation, namely at the collision site of heavy nuclei. 
After their collision, a quark-gluon plasma is believed to be formed ( see ref. [|J and refer- 
ences therein) with the temperature around 200 MeV. The plasma consists of mainly light 
quarks and antiquarks which may be considered to be ultrarelativistic, that is massless, 
particles under such a temperature. On the other hand, the charge of the original nuclei 
can not be destroyed in the process, but rather it surrounds the plasma putting it into a 
strong electric field. There is a possibility of light quarks running away in the very same 
way as the electrons in plasma. 

In this paper we consider this phenomenon in detail and show that due to the effects 
of relativistic mechanics the quarks will not run away and the only influence the electric 
field will have on them will be a slight distortion of the distribution function. 

To do it, we need to derive the kinetic equation of the quarks in quark-gluon plasma. 
The main interaction mechanism here is the exchange of a single gluon between quarks 
since at such high temperatures the perturbative QCD becomes applicable. The resulting 
interaction is very similar to the standard QED photon exchange, so the force between 
quarks will be of the same functional form as the electromagnetic force between electrons 
in plasma. This brings us to the conclusion that we should use the Landau collision integral 
(see 0) for the kinetic equation. 



In deriving the kinetic equation of the electrons in plasma, Landau noticed that the 
collisions with a small momentum transfer become important because of the slow decrease 
of the Coulomb force with distance. It allows to consider the collision process as a diffusion 
in the momentum space. 

Let us review the derivation of Landau Collision Integral following but not assuming 
that the plasma we are considering is nonrelativistic having in mind the application to 
quark-gluon plasma. 

The collision integral can be expressed in terms of the divergence of the flux of the 
particles of plasma 

4[ = _dsc L (1 , 

dt dp a 

where / is the distribution function, p a is the momentum and s a is the flux. Let us find 

that flux in terms of the distribution function. 

Let 

wf(p)f'(p')d 3 qd 3 p' 

be the number of collisions per unit time between the particles of the momentum p and p', 
q being the momentum transfer. We remember that we will be interested in collisions of 
high energy runaway particles with the thermal particles, so we denoted the distribution 
of high energy particles by / while denoting the Maxwell distribution by /'. 
We will write down w as a function of p, p', and q in the form 

n P+ 2' P ~2' q J- 
Here the first argument is half sum of the initial and final momenta of the first particle, 
while the second one is half sum of the initial and final momenta of the particle the first 
one collides with. 

The evident property of w is 

W ( P + f' p/ -f' q ) =M; ( P+ l' p/ -l'- q ) (2) 

Now we notice that by definition the flux s a (p) is the excess of particles which mo- 
mentum changes from some value smaller than p a to another one bigger than p a over 
the particles which momentum changes in the opposite direction, from bigger to smaller 
values. 



The total number of particles increasing their momentum is given by 

/ d\ J rfV r dp a w (p + |, p' - |, q) /(p)/'(p') 

while the total number of particles decreasing their momentum is 

/ d 3 q fd 3 p' I" dpaw(p+?, p'-?, -q)/(p + q)/V) 

Because of the relation (|2]) the functions w are the same in both expressions. We need 
to compute the difference of these two integrals. According to what was said above, the 
most important contribution to this expression comes from the small values of q. So we 
can expand everything in powers of q keeping only the first term of the expansion while 
the integral over p a can be replaced by multiplying by q a . We obtain 

f J* [ J* ' ( f „\ (ft xWPO fl, t\ df(p) \ . . 

/ d q / d pw(p,p,q) /(P)— 5-/ /(P) 9o9/3 (3) 

Jq a >0 J V ^ P/? / 

The integrand of @ is an even function of q due to (||) . So we can replace the integral 
over q a > by half the integral over all q. 

Now we can introduce the collision cross section by using 

cl 

wd 3 q = da , (4) 

where I is a well known invariant 



(P1P2) 2 -m\m\, 



Pi, P2 are the particle four momenta, m are their masses, and e are their energies. To 
make the calculations simpler we are going to write all the four- vectors and to do all the 
relativistic computations with the units where c = 1 while putting c explicitly in all the 
kinetic formulae. It should not result in any misunderstanding. 
After substituting ([|) into (|3|) we arrive at 

s a = j dp [f<P)-Qjr- ~ f (P )— J B aP (5) 

where 

1 f l 

B a /3 = o / 9a 9/3 cdo- (6) 

* J £1^2 



Thus far we were able to proceed without assuming anything about whether the 
particles are relativistic or not and what kind of interaction governs their collisions. Now 
our next task is to compute the quantities B which are of course dependent on what kind 
of plasma we deal with. From now on, we are going to consider the case of ultrarelativistic 
particles. 

First, we already saw that the distribution /' describes the bulk of quarks distributed 
according Maxwell distribution, while the distribution / is the one of the possible runaway 
quarks we will need to determine later. So 

/V)=Cexp(-|^) (7) 

where C is the normalization coefficient, c is the speed of light (the quarks are ultrarela- 
tivistic!) and T is the temperature. 

The tensorial structure of B a p can be determined by a very simple argument. If / 
were also a Maxwell distribution of the form ([7|) it would make the flux (|5|) to be zero. 
Let us substitute (0) into (|5|) for both f and /'. The integrand of (0) will then become 
proportional to 

l ff (?£ _ El) Bga 

T \p' p J 
On the other hand, it is better be zero. From here we see that B a p is a transverse tensor 
with respect to the difference of the unit vectors lying in the direction of the momenta of 
the colliding particles, 



B a /3 = ~B 5 a p —% (8) 



where B is just the trace of B a p, 

B = \U 2 —cda (9) 

2 J eie 2 

Another way of seeing the same is to notice that if the momentum transfer q is a very 
small vector in comparison with the momenta of colliding ultrarelativistic particles it will 
be orthogonal to the difference — — ^ due to the energy conservation and then (|8|) follows 
from (161). 
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Now let us proceed with computing B according to ([)]). The differential cross section 
of the electron-electron collision due to the photon exchange in the ultrarelativistic limit 
can be taken from, for example, and is equal to 

e 4 (3 + cos 2 #)) 2 , / N 

da = — ± ^—do (10) 

e 2 4sin 4 6> V ' 

where e is the energy of each of colliding particles in the center of mass reference frame and 
6 is the standard scattering angle in the same frame. There is an interference factor in ( TH]| ) 
which takes into account that the electrons are identical particles, but this factor is close 
to 1 at small scattering angles 0<1. The whole nature of Landau collision integral means 
only small scattering angles are important so the interference terms can be neglected. 

We will take the formula (|l0l) as the quark-quark collision cross section substituting 
a instead of e, a being the strong interaction effective coupling constant at the energies 
we consider (a ~ 0.1 . . .0.2). So 

a 2 (3 + cos 2 #)) 2 , / N 

da=^± ^^do 11 

e 2 4sin 4 6> V ; 



By taking the cross section in the form ([11]), we assume that the momentum of fast quarks 
is much larger than the Debye mass, meaning that the typical scattering is at distances 
much smaller than the Debye length. 

Now do is the spherical angle differential in the center of mass reference frame. But in 
(0) we integrate q 2 , which is the transferred momentum in the laboratory frame. In order 
to perform the integration, we have to express q 2 in terms of the angle 9. In general, that 
turns out to be a long algebraic formula which is difficult to work with. However, we can 
considerably simplify the calculations by using the fact that one of the colliding particles 
is a very fast quark, while the other is slow, even though both of them is ultrarelativistic 
(don't forget that we consider fast ranaway quarks in the "sea" of the slow Maxwell ones) . 
In this approximation it is rather easy to transform the ultrarelativistic momenta to the 
center of mass frame and then back to express the momentum transferred q in terms of 
the scattering angle. 
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Fig. 1: Quark scattering 



fig. 1 shows the 3-momenta of fast (P) and slow (p) quarks. We choose the frame 
of reference in such a way that the projections of P and p on the vertical axis Y are the 
same. Then the horizontal exis X will be the axis of the boost transformation we need to 
perform to pass to the center of mass reference frame. Since P ^> p, the angle between 
P and the X-axis is very small and the angle a between p and the axis X-axis almost 
coincides with the angle between P and p. If we denote p' to be the momentum of each 
quark in the center of mass frame, then from fig. 1 it follows that the projections p' x and 
p' y of p' to the axes X and Y are (here we choose the speed of light c = 1) 



Px 



P-Pv 



(12) 



Py—p sin a 
where v is the velocity of the boost. On the other hand, 

p cos a — pv 



-P'x 



Vl —v' 



So 



P + pcosa p . 

— — ~ 1 (1 — cos a) 

P + p P y ' 



(13) 



(14) 



Now let q' be the 3-momentum transfer in the center of mass frame. Its direction with 
respect to p' is defined by the scattering angles. The energy transfer in this frame is zero 
(the 3-momenta of the colliding particles just rotate without changing their magnitudes). 

q' 2 =4j9 /2 sin 2 - (15) 



9 being the scattering angle. We need to find the momentum transfer q in the laboratory 
frame while keeping it expressed in terms of the center of mass scattering angles. 

To do that, we compute the projections of the vector q' to the axes X, Y, and Z (Z 
is perpendicular to the plane of fig. 1) to obtain 

q' = p'(cos9 — l)(cos/3, sin (3, 0) + p' sin 9{— costp sin (3, cos ip cos /?, sin ip) (16) 

where (3 is the angle between the X-axis and p' while <p is the azimuthal scattering angle. 
One may check that ([16]) square gives exactly fll5p . 
The square of the four- vector q^ is clearly just 

2 il a il ■ 2 " 

q 4 = -q = -4p sin - 

Using that we transfer q' back to the laboratory frame to obtain 

12 2 

q 2 = -<& + « 2 = -<& + ' /j ' - 



» 2 « 2 (17) 

4p' 2 sin 2 — hp' 2 ((cos 6* — 1) cos /3 — sin 9 cos <z> sin j3) =■ 

2 1 — v 

Now we use that only the small scattering angles are important, that is 9 <S 1, to get 

1 — ir 
And the last step is to use that 



q 2 = p ' 2 9 2 [ 1 + cos 2 ip sin 2 /3- — ) (18) 



the formula (|TJ) and 



'2^/2 1-cos a 
P ~P X =P p o > (19) 



p' /2p 1 

sin/? ps — = sin a\ : (20) 

p x V P V 1 — cos a 

Substituting (|TJ), (0), and @IJ) to the (|TS|) we get 

> > 1 — cos a ( 9 sin 2 a \ , , 

- 2 - 2 pP ( 1 + cos 2 p- -2 (21) 

(1 — cos a) J 

This completes the task of expressing q in terms of scattering and azimuthal angles. 
Before proceeding further, we can average over the azimuthal angle to get 

., 1 — cos a / sin 2 a , 

0*P P n I 1 + TZ 72 I ( 22 ) 

2(1 — cos a) 



Now we can plug everything into the (|J). Keeping in mind that the most important 
angles are the small ones, ^Cl, and using 

/ 





eif2 


1 — cos a 
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i sin 2 a \ 
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, 2(1 — cos a) I 


J o 



we arrive at 

B = ^ (1 -°" o) i 1 + i(ir^jyf (23) 

The integral over 9 turns out to be logarithmically divergent. The same thing happens 
to the nonrelativistic plasma as well so we should not be discouraged. In fact, it should be 
cut off at some small angles due to the finite density of plasma and so, the impossibility 
of very small angle scattering, and it large angles because (23) is applicable only when 
<C 1. Within the logarithmic accuracy, it is enough to replace the divergent integral by 
some number L called the Coulomb logarithm. Then B becomes 

„ 47ra 2 L ,„ . / sin 2 a \ ,„„. 

B = (l-cosa)lH o 24 

c J \ 2(l-cosa) 2 y 

B has a simple meaning of the average momentum square transfer in a collision. We 
see that B does not depend on the momenta of the colliding particles at all unlike the 
nonrelativistic case where B decreases with the momentum as 1/P (see ref. ||). As we 
will see later, that changes completely the behavior of fast quarks in comparison with fast 
electrons. While fast electrons run away in the direction of the electric field, fast quarks 
do not. 

Now we are in the position to derive the final expression for the flux from @. To 
do that, we substitute (|24l ) with (||) into (|f|). Then we remember that the function f'(p') 
is a Maxwell distribution (|7|) while the function f(p) is the unknown distribution of the 
ranaway quarks. And we integrate out the momentum p' of the equation (^). There are 
two terms in (H) we have to compute. Let us begin with the first one. 



J d p ^pV a " f (P } 



We introduce the unit vectors n = p/p and n' = p'/p'. In terms of them the expression 
becomes 



2 LN df I do' /3 1 , 

-Zna — — / - — n^n~ 

c dp p J 4n \2 2 7 7 



(n a -n' a )(n0-n'g) 



2n^n' y 



where the integration ^- gives the average over the spherical angle of n' and we used that 



the integration over the magnitute of p' gives us just the density of quark plasma: 

N = 4tv I p' 2 dp'f'(p') 



To compute the integral over the spherical angle we notice that the only vector the 
result may depend on is n. So 



do' . 

An V 2 ~ 2 nin i 



lap 



(ng-n'Jjnp-n'p) 
2 - 2n 7 n; 



= CxSa/3 + C 2 n a ni3 



By taking the trace of this expression and by multiplying it by np we immediately obtain 



J Jn~ \2 2 
Therefore we are left with 



4^l2~2 n ^ 



lap 



(n a -n' a )(np-n' p ) 
2 - 2n 7 n 7 



7 A 1 

^o a p - -n a np 



2 LN df (1 1 

-2-kol — -b a/ 3 - -n a np 

c opp \6 2 

We will see later that we need just the radial component of that expression. It 

4na 2 LN df 



is 



(25) 



3c dp 

We don't have to compute the second term of (^) which is proportional to /. The expression 
for the radial component of the flux should look like 



C^M + B/(P) 



(26) 



It turns out we do not need to compute both C and D in ( p6| ) because they are related. 
The Maxwell distribution must make the flux (|26|) be zero. That is 



C 



or 






(27) 



So it is enough to compute the coefficient in front of the term of (|5|) with the derivative of 
/. But it had been already computed and given by (|25|). Finally we obtain 

4a 2 LN fdf{p) c 
3c V dp T' 



-• = : '" I — — + ™/(p) 



(28) 
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The main difference between the flux of ultrarelativistic particles we have just obtained 
with that of classical particles found in |5| is the momentum dependence of the "classic" 
flux. It is Selassie oc --gr while relativistic flux does not depend on the momentum at all. 
That strong collision flux will eventually prevent particles from accelerating in the direction 
of the electric field. 

Now, with the expression for the flux in our hands, we can proceed to discuss the 
kinetic equation. 

The general form of the kinetic equation is easy to write down as soon as the expression 
for the flux is known 

| +eE | + diVpS = (29) 

where e is the charge of the particular quark and E is the magnitute of the electric field. 

In general the function / depends on both the coordinates and the momenta. The 
geometry of our problem is spherical. We have a sphere of quark-gluon plasma surrounded 
by the electric charge which means spatial nonuniformity of the electric field. Nevertheless, 
the basic question whether the quarks are going to run away or not does not depend on 
any geometry. Since it is much easier to consider the quarks in the homogeneous electric 
field, we are going to solve that problem. Our spherical problem cannot give rise to a 
totally different behavior, qualitatively the answer will be the same. 

We conjecture now that there is such a phenomenon as runaway quarks. That means 
the solution of the kinetic equation will give us the function /(p) which will describe the 
runaway quarks moving in the direction of the electric field. To find the flux of such 



particles we can average the kinetic equation over the directions, following ||. ( |29|) gives 



us 



% + eE 4 r P 2 f + 4rP 2 ^ = o (30) 

at p z op p z op 



while the stationary solution of kinetic equation satisfies, with the aid of (BS| 

"4a 



,-, d 9 „ d 9 
eE^—p z f - -yttV 
p z op p z op 



x 2 LN fdf{p) c \ 



(31) 



This is the fundamental equation we have been deriving. 

Let us try to find the constant flux stationary solutions. They are given by 



i- 



J 3c V dp T JyF> 
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n (32) 



where n is the flux. We can easily check that if the flux n is not equal to zero, the solution 
of (^2j) will always be divergent at p — > making it impossible to find such a solution which 
describes a constant flux at large p while smoothly connected with Maxwell distribution 
at small values of p. That means there are no runaway quarks, as we have promised. 
We can try to analyze the zero flux solution, 

f = e -{T- 47ra 2 LN )p_ (33) 

It gives us a slightly distorted Maxwell distribution. 

Let us estimate the terms in ( |33D to see how small the correction is. T is of the order 
200 Mev. E can be estimated as being created by a charge about 100 x e at the distance 
of about oc 10 _13 cm. The density iV is the number of quarks per unit volume. Let us also 
estimate the number of quark as 10 2 then iV oc 10 41 cm -3 . That gives us 

4ocl0 14 S 



T g cm 

and 

3eEc „„•,„ s 



oclO 
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ATia 2 N g cm 

So it is difficult to say if the electric field is strong enough to produce a noticible 
distortion of the distribution function on the basis of this very rough estimation. Let us, 
however, notice that the whole derivation was based on the idea that only fast quarks are 
distorted while the distribution ( |33| ) would predict an overall distortion. So most probably 
( P^D is unphysical. 

In conclusion we would like to say that we have shown the phenomenon of the run- 
away quarks does not exist due to the effects of the relativistic kinematics. One can also 
undertake a similar calculations with c-quarks which are massive enough to be almost 
nonrelativistic at our temperatures. We will have the problem of nonrelativistic c quarks 
trying to run away in the sea of the ultrarelativistic u and <i-quarks. This calculation 
turns out to be similar to the one given in this paper. We will not give the details of that 
calculation here, but we will remark that the result is the same, the heavy quarks do not 
"run away" as well. 



The author is grateful to I. Kogan for posing the problem and for many discussions on 
its subject and also to the Physics Department of Melbourne University where this work 
had been completed for their hospitality. 
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